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1 Quantum manifestations of the classical 
stochasticity 

1.1 The quantum chaos problem 

Essential progress in the understanding of the nonlinear dynamics of classical systems 
stimulated numerous attempts to include the conception of the stochasticity in quantum 
mechanics [0, [|, [H], [TTJ, F2, F3, [TJJ. The essence of the problem consists in the fact that 
the energy spectrum of any quantum system, performing the finite motion, is discrete 
and therefore its evolution is quasiperiodic. At the same time, the correspondence prin- 
ciple requires the possibility of limit transition to classical mechanics, containing not 
only regular solutions but also chaotic ones. The total solving of this problem is too 
complicated one, therefore its restricted variant is of considerable interest. We are inter- 
ested in the search of peculiarities of behaviour of quantum systems the classical analogs 
of which reveal chaotic behaviour. Such peculiarities are usually called the quantum 
manifestations of the classical stochasticity (QMCS). 

The following objects can be used as the objects of search of the QMCS for an 
autonomous systems 
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1) the energy spectra 

2) the stationary wave functions 

3) the wave packets. 

A priori the manifestations of the QMCS can be expected both in the form of some 
peculiarities of concrete stationary state and in the whole group of close states in energy. 
Of course, it is not excepted the possibility that such alternative does not exist at all, 
i.e. the traces of the classical chaos can be detected both in the properties of separate 
states and so in their sets. 

In this part of the review we will represent the results relating to the QMCS in the 
dynamics of quadrupole oscillations. 

1.2 Numerical procedure 

The quantum scaled Hamiltonian describing the quadrupole surface oscillations has the 
following form 

H=^(A x + A v ) + U i (x,y,W) (1) 

where the scaled Planck's constant is H = fyt , h\ = me ^oj > an d scaling parameters 
hoijoi and deformation potential Ui(x,y,W) are determined by the expressions (21.1) 
and (23.1). We are going to study the peculiarities of structure of energy spectra and 
wave functions in each of the following intervals 

the first regular region Ri : < E < E cr \ 

the chaotic region C : E crl < E < E cr2 

the second regular region R 2 '■ E > E cr2 . 

The critical energies E cr i and E cr2 were determined above. The main numerical 
calculations were performed for one- well potentials (0 < W < 16 ). In this case the 
critical energies are determined by the relations (60.1), (61.1). 

At fixed topology of potential surface (W = const) , the scaled Planck's constant hi 
is the unique free parameter in the quantum Hamiltonian (fj). In the study of the con- 
crete region of energies, corresponding to the certain type of classical motion (R 1 , C, R 2 ) 
, the choice of hi is dictated by the possibility of attainment the required degree of valid- 
ity of semiclassical approximation (large quantum numbers) conserving the precision of 
calculations of spectrum and wave functions (the limitation of the possibility of diago- 
nalization of large dimension matrices). For the original nonreduced Hamiltonian (16.1) 
the variation of Hi is equivalent to the variation of the critical energies E cr \ and E cr2 (i.e. 
parameters a,b,c ) or to the adequate choice of Planck's constant H . 
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Table 3.2.1 





A 1 


A 2 






67 


70 


48 


m 


2346 


2556 


1225 



The procedure of diagonalization was used for the determination of energy spectrum 
and eigenf unctions of the Hamiltonian ([]]). As the basis we chose the simple combinations 
of the eigenfunctions of the two-dimensional harmonic oscillator with equal frequencies 
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\NLj) = ^ (\NL) + j \N, -L)) , j = ±1 
iV = 0,1,2...; L = N,N — 2...1 or 
normalized according to the condition 

(NLj | N'L'f) = l=2 s ^5 jf 5 NN/ 8 LL , (3) 

and P LJ = j Mod ( L & . 

The symmetry of the considered Hamiltonian (C^ v ) leads to the block structure of 
matrix (N'L'j'\H\NLj) , which consists of four independent submatrices. It allows 
to carry out diagonalization of each submatrix separately, obtaining in this case four 
independent sets of states, which according to the classification used in the theory of 
groups are the following states 

A\ - type [mode(L, 3) = 0, j = 1 , including L = ] 

A 2 - type [mode(L, 3) —0,j — —1, L^O] 

and twice degenerated states 

E - type [ mode (L, 3) ^ 0, j = ±1 ] 

The possibility of separate diagonalization of submatrices of definite type allows 
essentially to increase available dimension of basis in numerical calculations (later for 
brief, we shall call the dimension of basis as the dimension n of separate submatrix). 
Table 3.2.1 shows the maximum oscillator numbersAo of the basis and the dimension 
m of total matrix (N'L'j'\ H \NLj) for the above enumerated types of submatrices with 
the dimension n = 408 
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Figure 18: Energy spectra of Hamiltonian (jl|) as function of the frequency ujq of the 
oscillator basis (W = 13) . The upper and lower parts of this figure show this dependence 
for the basis dimensions n = 198 and n = 408 , respectively. 

Expansion of wave eigenf unctions in basis (0) in polar coordinates r, ip 

(r, V >\E k ) = J2C^l(r, V \NLj) (4) 

NL 

where 

(5) 

(V^r) L e- r2 /2M {-^A L + l, ^ } 

Here M{ } is degenerated hypergeometric function and u> - frequency of the oscillator 
basis. The finite dimension of the basis used in calculations leads to the dependence of the 
position of level from the oscillator frequency c^o • Such dependence for the dimensions 
of basis n = 198 and n = 408 (W = 13 ) is represented in Fig. 18. 

The important moment of numerical calculations is the choice of the optimal fre- 
quency Uon opt , determined from the condition that the energies of levels, situated in the 



(r,cp\ NL) 



■N 



-iLtp j 

27 L\ 



(¥)' 



5 



interval of our interest, have minimum value. Such optimization is essentially impor- 
tant for the regions C and R2 , where the position of the levels more depends on the 
frequency of basis. The optimal frequency depends both on dimension of the basis and 
on the length of investigated energy interval. As is shown in Fig. 18 it is possible to 
choose the unique frequency ujQ Pt for intervals including decades or more levels just at 
the dimension of the basis n = 408 . These levels are selected in Fig. 18 by dotted lines. 
The large value of u)^ for the region R 2 (u)° ] pt ~ 8) is caused by the fact that at energies 
E > E cr 2 the deformation potential essentially differs from harmonic oscillator potential. 

The dimension of the basis used in our calculations was determined for the reason of 
acceptability of the calculation time and the calculation precision of the position of levels 
in considered region. The results of investigation of saturation in basis are represented 
in Fig. 19, i.e. dependence of the results of spectrum calculation from the dimension 
of basis. As it is shown in Fig. 19 for the regions R\ and R2 acceptable precision of 
calculations of levels with ordinal numbers 100 — 200 is reached at the dimension of the 
basis n ~ 400 , while for the chaotic region C dimension of the basis required for the 
attainment of the same precision considerably larger (n ~ 700) . It is a well known result 
|3]. |(|: the value of saturation in basis is sensitive to the type of classical motion. 

Matrix diagonalization is attractive for treating Hamiltonians that do not differ 
greatly from Hamiltonians with known eigenf unctions. This numerical procedure be- 
comes less attractive (or even not effective at all) at the transition to the PES of com- 
plicated topology, when required Hamiltonian matrix is of high order. In this case the 
alternative to the diagonalization may become the so-called spectral method |16|] , which 



utilizes numerical solutions to the time-dependent Schrodinger equation. The spectral 
method was developed earlier for determining the eigenvalues and eigenf unctions for the 



modes of optical waveguides from numerical solutions of the paraxial wave equation fT7 



Feit, Fleck and Steiger |I6| were able to apply the previously developed methodology to 
quantum mechanical problems with little change, as the latter equation is identical to 
the Schrodinger equation. 

The spectral method requires computation of the correlation function 

P 1 (t) = (if>(r,0)\if)(r,t)) (6) 

where if> {r, t) represents a numerical solution of the time- dependent Schrodinger equa- 
tion, and if) (r, 0) is the wave function at t = . The solution if) (r, t) can be accurately 
generated with the help of the split operator FFT (Fast Fourier Transformation) method 



[ i~8| . The numerical FFT of P\(t) , or P\{E) , displays a set of sharp local maxima for 
E = E n , where E n are the desired energy eigenvalues. With the aid of lineshape fit- 
ting techniques, both the positions and heights of these resonances can be determined 
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Figure 19: Energy levels as function of the basis dimension n . The numbers on the 
curves are the level numbers k . 
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with high accuracy. The former yield the eigenvalues and the latter the weights of the 
stationary states that compose the wave packet once the eigenvalues are known, the 
corresponding eigenfunctions can be computed by numerically evaluating the integrals 

T 

t/j (f, E n ) = J^j (f, t) W (t) exp (iE n t) dt (7) 
o 

where T is the time encompassed by the calculation, and W(t) is a window function, 

W(t) = l- cos 27Tt / T if < t < T 

(8) 

= if t>T w 

Since the spectral method is fundamentally based on numerical solutions of a time- 
dependent differential equation, its implementation is always straightforward. No special 
ad hoc selection of basis function is required, nor is it necessary for the potential to have 
a special analytic form. The spectral method is in principle applicable to problems 
involving any number of dimensions. 



1.3 Quantization by the normal form 

In this section we calculate a semiclassical approximation to an energy spectrum of the 
Hamiltonian of qudrupole oscillations (0) |19| (here it will be more suitable for us to 
use the nonscale version of the Hamiltonian (16.1) with m — 1 ) by quantization the 
normal form and compare results to exact quantum mechanical calculations. By the 
exact spectrum we mean the spectrum obtained by the direct numerical calculations, 
for example, by the diagonalization of the Hamiltonian on the reasonably chosen basis. 
Quantization the incommensurable case is straightforward. Since the normal form can 
be expressed entirely in terms of action variables, we merely have to replace those action 

variables by an appropriate multiple of h, I v — > h (n u + l/gj and the result is a power 

series in the quantum numbers, 

E nin2 =T(n 1 + 1 / 2 ,n 2 + 1 / 2 ) (9) 

The commensurable case is less straightforward. Commensurability of frequencies leads 
to the appearance of angle variables in the normal form 

r = r(/ 1 ,/ 2>¥>2 ) (io) 
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As before, the action I\ we can replace on h \ n± + I/2J . In this case the initial two- 
dimensional problem is reduced to one-dimensional problem, the quantization of which 
can be performed with the help of WKB method. 

The procedure of quantization by normal form we shall begin with canonical trans- 
formation 

q\ = \ (~Qi + Q2 + P1- P2) 
Q2 = l (Qi + Q2 + P1 + P2) 

\ (11) 

Pi = 2 (Gi -Q2 + P1- P2) 
V2 = \ (Qi -Q2-P1- Pi) 

where variables (p, q) provide reduction of the harmonic part of the initial Hamiltonian 
(16.1) to the form (77.1). In the new variables Q(Qi, Q2) and P{Pi, P2) the Hamiltonian 
of quadrupole oscillations is the following 



K (Q, p) = (Q, P)+J2 KU) (* P) ( 12 ) 

i>2 

where 

K( 2 )(Q,P) =*(QiPi + Q 2 P 2 ) (13) 

and K&> are homogeneous polynomials in variables Qi and Pj of degree j. Each member 
K^) of the Hamiltonian (0) is reduced to the normal form according to the procedure 
describing in the section 2.5(part 1). However this procedure is sufficiently simplified 
due to the diagonal form of the operator D in the variables (P, Q) . The classical normal 
form is the sum of the polynomials of special form in Qi and Pj . For obtaining its 
quantum analog we can use the Weyl's heuristic rule of correspondence [20, 21 

1 n n\ 

pn Q m = Q mpn ^ _ £ _ __^ Qm pn-l {u) 



1=0 



The operators p and Qi are determined by formulae (|TT|) , in which by pi and qi we 
should mean operators of impulse and coordinate with usual rules of commutation, from 
which it follows that 



Pfe> Qi 



MM = 1,2) (15) 
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The operators Pj , Qi allow to introduce the full orthonormalized basis 



\NL) 



N + L 



N -L 



2 A N-L A iV+L 

Q 2 2 Qi 2 |o> 



where vacuum state |0) is determined by 

A |o) = A |o) = 



(16) 



(17) 



The principal quantum number N = 0,1, ... and angular momentum number L at given 
N is equal ±iV, ±(iV — 2)...0 or 1. The action of the operators Qi and Pj on the basis 
\NL) is 



Qi |iVL) 
Q 2 |iVL) 



iV + L + 2 



|JV+l,L + l),Pi|i\TL) 



iV-l,L-l) (18) 



N-L + 2 



N + 1,L-1), P 2 \NL) 



N-L 



N-1,L + 1) 



Note that the constructed basis is oscillator one 

K (2) \NL) = (QiA + g 2 A + l) \NL) = (N + 1) |iVL) (19) 
The classical normal form for considering Hamiltonian is (through fourth degree) 

r = r {2) + r (4) = 

i [{Q1P1 + Q2P2) + 6 7e ( p iQl + p iQl - i2QiQ 2 PiP 2 ) + 

+c (PlQ\ + P^Ql + 4g 1 g 2 p 1 p 2 )] (20) 

The quantum normal form T reconstructs from the classical one with the help of Weyl's 
rule ( |14D and Dirac's rule of correspondence T — > hT 



QxPi + Q2P2 + 1 



b 2 

+ 6 



QiPi) 2 + (Q2P2) 2 -5(Q 1 P 1 + Q 2 P 2 ) ~12Q 1 P 1 Q 2 P 2 -2 



+c 



QxPi ) + (Q2P2) 2 + 3 ( QxPx + g 2 p 2 ) + 4g 1 p 1 g 2 p 2 + 2 



(21) 
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It is easy to see, that basis vectors \NL) are eigenvectors for quantum normal form 
Therefore we get the simple analytical formula for the energy spectrum in the fourth 
approximation. 

E (N, L) = N + 1 + ^ [lL 2 - 5 (N + l) 2 + l] + | [3 (N + l) 2 - L 2 + l] (22) 



Assuming c = in the formula ( [22]) we get approximate spectrum of the Henon-Heiles 
Hamiltonian, which up to the constant shift coincides with the spectra obtained by other 
methods. The last property, apparently, is connected with the ambiguity of quantization 
method of the Hamiltonian. Each level in the energy spectrum (22) is double degenerated 
with respect to sign of angular momentum number L , whereas for the levels of the exact 
Hamiltonian with L = 3k(k = 1,2, ...) the degeneration must be taken off. Inclusion in 
the normal form the members of higher degree leads to the taking off the degeneration. 
Really 

T (6) = + y b 2 c + 2c 2 ) (P?Ql + PiQl) + 

+ (^& 4 - y b 2 c + 15c 2 ) P l P 2 Q l Q 2 (PxQx + P 2 Q 2 ) + (23) 

The last term in this relation connects the states with equal and the states with 
V — L = ±6 . However in this approximation the basis vectors \NL) ([16]) will be no 
longer the eigenvectors of the normal form. Therefore for the calculation of the energy 
spectrum it is required either additional diagonalization of the engaging states or the 



calculations by the theory of perturbation [20]. How well does the energy spectrum 
fl2"2|) , obtained by quantization of classical normal form in the fourth approximation, 
reproduces exact quantum spectrum of the Hamiltonian of quadrupole oscillations? In 
the Table 3.3.1 approximate spectrum E(NL) is compared with the spectrum E ex , 
obtained by diagonalization on oscillator basis. 

The dimension of submatrices of type A ly A 2 ,E are respectively equal 560,560 and 
415 and provide sufficient precision of calculations for the first 250 — 300 levels. The 
parameters of the Hamiltonian (16.1) are chosen so that E cr \ = 90 . As it can be seen 
from the Table 3.3.1 in region, where classical motion is regular, the quantum normal 
form reproduces the energy spectrum with the relative error AE/E ~ 0.01% . It seems 
naturally to expect that in the neighbourhood of critical energy of the transition to chaos, 
where the destruction of the approximate integrals of motion takes place and with the 
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Table 3.3.1. Energy levels of Hamiltonians (16.1) (6 = 0.04416, c = 0.00015, E cr = 90) 
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Figure 20: Difference between exact quantum energy levels and quasiclassical levels (]22f) . 
The line is the mean level spacing, (a) for Henon-Heiles Hamiltonian with parameter 
b = 0.1 (E cr = 8.3; E$ = 16.7) ; (b) for Hamiltonian of quadrupole oscillations with 
parameters: W = 13; a = 1; c = 0.00135 (E cr = 10). 
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help of which the semi classical spectrum is built, the agreement of the last one with 
the exact spectrum must essentially be getting worse. Let us analyze this effect on the 
example of the spectrum of quadrupole oscillations. The difference AE = \E (M) — E ex \ 
as the function of energy calculated for the case W = 13, E cr = 10 is represented in 
Fig. 20(b). We can see that in the region of energies, where the classical motion is 



regular, the approximate spectrum (22) reproduces the exact one rather well. At the 
transition to chaotic region the difference increases sharply. The similar situation takes 
place and for the Henon-Heiles Hamiltonian (Fig. 20(a)). 

1.4 R — C — R transition and statistical properties of energy 
spectrum of the QON 

Important correlations between peculiarities of classical dynamics and structure of quan- 
tum energy spectrum can be obtained in investigation of statistical properties of level 
sequences. We shall be interested in local properties of spectrum, i.e. in deviation in 
distribution of levels from mean values and in fluctuations. Why should we address to 
local characteristics of spectrum? The matter is, that the global characteristics like the 
numbers of states N(E) or the smoothed density of levels p{E) are too rough. At the 
same time, such local characteristics, as the function of the nearest-neighbour spacing 
distribution between levels is very sensitive to the properties of potential and to the 
shape of boundary. It is sufficient, for example, to bend slightly one of the walls in 
square to make it scattering, then classical trajectories in such system become chaotic. 
In this case the essential rearrangement of functions of the nearest-neighbour spacing 
distribution between levels happens, though the number of states N(E) changes slightly 
or doesn't change at all. 



Relying on rather simple reasons [p2| , let's try to construct one of the local character- 
istics of quantum spectrum - the function of the nearest-neighbour spacing distribution 
between levels p(s) . For a random sequence the probability that a level will be in the 
small interval (E + s, E + s + ds) , proportional, of course, to ds , will be independent 
of whether or not there is a level at E . This will be modified if we introduce level 
interaction. Given a level at E , let the probability that the next level) (s ^ 0) be in 
(E+s, E + s + ds) be P(s)ds . Then for P(s) , the nearest-neighbour spacing distribution 
(NNSD), we have 

P{S)dS = P{ledS\0eS)P{0eS) (24) 

where P(n G S) is the probability that the interval of length S contains n levels and the 
P(n G dS\m G S) is conditional probability that the interval of length dS contains n 
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levels, when that of length S containes m evels. The second factor in ( P4|) is J P (x) dx 

s 

, the probability that the spacing is larger than S , while the first one will be dS times 
a function of S,r w (S) , depending explicitly on the choice, 1 and , of the discrete 
variables n, m . Then 



P(S) =r 10 (S) J P(x)dx (25) 
s 

which we can solve easily to find 

P (S) = cr 10 exp \ - I no (x) dx 1 (26) 



The Poisson law follows if we take rio(S) = 1/D (the absence of correlations between 
level positions), where D is the mean local spacing so that 1/D is the density of levels. 
Wigner's law follows from the assumption of a linear repulsion, defined by r w (S) = aS 
. The arbitrary constants are determined by conditions 



P (x) dx = 1, J xP (x) dx = D (27) 
Finally, we find for the Poisson and Wigner cases, respectively 

P(S) = l/Dexp(- S /D), S^O (28) 



P ^ = 2^ exp S>0 (29) 



The second distribution displays the repulsion explicitly since P (0) = , in contrast to 
the Poisson form, which has a maximum at S = (level clusterization). 

To get a first idea about origin of the level repulsion let us consider |22[ the Hamil- 
tonian as defined with respect to some fixed basis by its matrix elements. The repulsion 
may be regarded as arising from the fact that the subspace for which the corresponding 
spectrum has a degeneracy is of a dimensionality less by two than that of the general 
matrix-element space, so that in some sense a degeneracy is "unlikely". Alternatively 



[23], if we think of the matrix elements as functions of a parameter a , we cannot, in 
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general, force a crossing by varying a but we must instead take the matrix elements 
as functions of at least two parameters which are independently varied. In the one- 
parameter case one will find that, if two levels approach each other as a is varied, then 
instead of crossing they will turn away as if repelled. 

There is a principal difficulty with the derivation of (|29|) . Why should we assume a 
linear repulsion? Although there are some plausibility arguments for this form, the result 
cannot be correct for every system. Furthermore only a probability arguments cannot 
explain the nature of level repulsion. However the situation changes, if we intend that 
statistical properties of the sequence of levels of real physical system are equivalent to 
the sequences of eigenvalues of the ensemble of random matrices of a definite symmetry. 



The theory, based on this hypothesis, was completed in the sixties piU , The final result 
for the function NNSD is 

P {S) « S a exp (-PS 2 ) (30) 

The critical index a determining the behaviour of the distribution function at S —>■ 
depends on the symmetry of the ensemble of matrices. This symmetry is determined 
by the properties of the physical system, the spectrum statistics of which we want to 
reproduce. If the system is invariant relative to time reversion, then the corresponding 
ensemble is Gaussian orthogonal ensemble (GOA). For the system assuming the violation 
of the invariance relative to time reversion the Gaussian unitary ensemble is associated. 
Finally, the symplectic ensemble of random matrices corresponds to the Hamiltonian 
of more complicated structure H = H + hr , where = , h K = h* K = — , 
(K = 1,2,3; <jk - Pauli matrices). The critical index a in the formula ( PU| ) is equal 
to: a = 1 for the GOA, a = 2 for the unitary ensembles and a = 4 for the symplectic 
ensembles. 

The predictions of the statistical theory of levels (especially for the GOA) were com- 
pared in details with all available sets of nuclear data [^, |24|. No essential deviations 
between the theory and data base have been revealed. Similar comparisons have been 
realized for atomic spectra. And here, a good agreement has been revealed with the 
predictions of the GOA, though the number of processed data was essentially less, than 
in the nuclear spectroscopy. 

If for the complicated systems (atomic nucleus, many-electron atom) we can give seri- 
ous arguments in favour of the hypothesis of the equivalence of the statistical properties 
of the spectrum and the sequence eigenvalues of the ensemble of the random matrices, so 
its generalization does not seem natural in the case of the systems with a small number 
of degrees of freedom. 

A radically new and sufficiently universal approach to the problem of the statistical 
properties of energy spectra may be developed on the basis of a nonlinear theory of dy- 
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namic systems. The numerical calculations ^ ^3], ^ |57|. supported by sophisticated 



theoretical considerations |l|, |3|, £|, ||, [XT], [13], |29], |30|, ^TJ , show that the important universal 



peculiarity of the energy spectra of the systems, that are chaotic in the classical limit, 
is the phenomenon of the level repulsion; while the systems, whose dynamics is regular 
in the classical limit, are characterized by the level clusterization. This statement is 



sometimes called [25] as the hypothesis of the universal character of the fluctuations of 
energy spectra. 

Among the systems, which spectra were subjected to detail numerical analysis, the 
central place is occupied by two-dimensional billiards (free particle moving on the plane 
inside of some region and subjecting to elastic reflection on the boundary). One of 
two extremal situations can be realized for the billiards with the definite shape of the 
boundaries: integrable or nonintegrable. The angular momentum is the second integral 
(except energy) in the circular billiard and such a system is integrable. The billiard like 
"stadium" is one of the simplest stochastic systems [p7|]. 

The NNSD for the integrable system (the circular billiard) is well approximated 
by Poisson distribution and the variance is the linear function of the considered energy 
interval, that is in complete correspondence with the hypothesis of the universal character 
of the fluctuations of energy spectra. In the nonintegrable case ("stadium"), the effect 
repulsion of levels and slow growth of the variance, caused by the rigidity of corresponding 
spectrum, are observed. 

Measure of rigidity is the statistic A 3 of Dyson and Mehta 

x+L 

A 3 (L; x) = ^rMin A , B J [n (e) - As - B} 2 dx (31) 

X 

which determines the least-square deviation of the staircase representing the cumulative 
density n (e) from the best straight line fitting it in any interval [x, x + L] . The most 
perfectly rigid spectrum is the picket fence with all spacing equal (for instance, the 
one-dimensional harmonic oscillator spectrum), therefore maximally correlated, for which 
A 3 (L) = 1/12 , whereas, at the opposite, the Poisson spectrum has a very large average 
value of A 3 (A 3 = L/15) , reflecting strong fluctuations around the mean level density. 

In contrast to billiards, where the character of motion does not depend on energy, 
the Hamiltonian systems of general position are the systems with the separable phase 
space, which contain both the regions, where the motion is stochastic, and the islands 
of stability. How is this circumstance reflected in statistical properties of spectrum? 
Berry and Robnik |3(J and independently Bogomolny J33|, basing on the semiclassical 



arguments, showed that NNSD for such system represents the independent superposition 
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Table 3.4.1 





Ri 


C 


R2 


W = 13 


E < 8 • 10" 5 


8 • 10" 5 < E < 8.4 • 10" 2 


E > 8.4- 1(T 2 


h 


h= 3.2 • 1(T 6 


ft = 1.6 • 1(T 4 


h = i.i- lo- 1 



of the Poisson distribution with the relative weight \i , determining by the part of the 
phase space with regular motion, and the Wigner distribution with the relative weight 
//(// + p, = 1) , determining by the part of the phase space with chaotic motion 

P (x) = jJ 2 exp (—fix) erfc (^-^-Jix^j + ^2//// + ^ffix^j exp ^— \ix — — p, 2 x 2 j (32) 



The expression (|52]) represents the interpolated formula between the Poisson ( PEQ and 
Wigner (p9|) distributions. 

Now let us go to the statistical properties of the spectrum of the Hamiltonian of the 
quadrupole oscillations (p. We intend to study the evolution of these properties in the 
process of R—C—R transition, i.e. in each of energetic intervals Ri, C, R2 . In this study 
we are restricted by the case of one- well potentials ( W < 16 ). At the fixed topology of 
potential surface (W = const ) the unique free parameter of quantum Hamiltonian (JJ) 
is the scaled Planck's constant hi . In the study of the concrete energetic interval (R\, C 
or i?2 ), corresponding to definite type of classical motion, the choice of % is dictated by 
the possibility of attainment of the necessary statistical assurance (the required number 
of levels in investigated interval) with conservation of precision of spectrum calculation 
(restrictions to possibility of diagonalization of matrices of large dimension). Let us 
notice, that for the initial nonreduced Hamiltonian (16.1) the variation of hi is equivalent 
to the variation of the parameters a; b; c (i.e. variation of the critical energies E cr \ and 
E cr 2 ) or to the adequate choice of initial Planck's constant h\ = h . The values of 
the scaled Planck's constant h% = h , represented in Table 3.4.1, allow us to obtain in 
corresponding energy intervals some hundreds of levels with the precision better than 
1%. 



Two additional comments are necessary before proceeding to results. We first recall 
that, to get rid of spurious effects of the local properties due to variation of the density, 
one has to work at constant density on the average. For this purpose one can 'unfold' 
the original spectrum [p2]| ) i.e. map the spectrum of eigenvalues {Ei} onto the spectrum 
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{ei} through 



Si = N (Et) (33) 
Here the smoothed cumulative density N (E) , 

E 

N(E)= I p (E') dE' (34) 



and p (E) is the smoothed level density. In what follows we will take as energy unit the 
average spacing x between two adjacent levels of the unfolded spectrum 

x = Si = ({e i+1 - £i)) (35) 

Meaning of unfolding is understood by the following example []34[| . There are a few 
spacings between states of the same(J 7r ,T) in the ground state region known for many 
nuclei. Together they constitute a large sample of experimental spacing data. The only 
trouble is that the underlying scaling parameter D varies from nucleus to nucleus. As 
a result, we cannot make statistical studies unless we have a model for the variation of 
D as a function of nucleon number A . This is a kind of unfolding except, instead of 
energy, we must remove the A dependence of D . 

In the second place, we remark that our above discussion about Poisson of Wigner 
level statistics applies only for a pure sequence, i.e., one whose levels all have the same 
values of the exact quantum number. More precisely, a pure sequence represents a 
set of levels relating to one and the same nonreducible representation of the group of 



symmetry of considered Hamiltonian. In the case of mixed sequences |22j the level 
repulsion is moderated by the vanishing of Hamiltonian matrix elements connecting two 
different symmetry states, the spectra have moved towards random, and NNSD towards 
Poisson. The results of investigation of correlations between statistical properties of 
energy spectrum of the Hamiltonian QON and the character of classical motion are 
represented in Fig. 21 J55|| . The typical Poincare surfaces of section, showing the change 
of classical regimes of motion, are represented at the left of Figure, and the statistical 
characteristics of spectrum (logarithm NNSD and dispersion), illustrating the evolution 
of the latter in the process of the R — C — R transition, are represented at the right. 
At the construction of statistical characteristics, a purity of sequence is provided by 
using only those levels, which are relative to definite nonreducible representation of C^ v 
group (the levels of E -type were used for results represented in Fig. 21; the statistical 
characteristics of levels of A\ and A2 - have similar form). 
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Figure 21: Correlation between the character of classical motion and statistical properties 
of energy spectra in R — C — R transition for the Hamiltonian of quadrupole oscillations. 
On the left: - Poincare sections, in the middle - distribution function of the nearest- 
neighbour spacings, on the right - the variance. Below - the first regular range (Ri) , in 
the middle - the chaotic range, above - the second regular range (R 2 ) ■ 
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Both the NNSD and the variance well correspond to the predictions of GOE for 
the chaotic (C) region . The logarithmic scale for NNSD p(x) is suitable to trace this 
correspondence at large x . For regular regions R\ and R2 the distribution function, 
in the same scale, according to the hypothesis of the universal character of fluctuations 
of energy spectra, must be represented by the straight line (the logarithm of Poisson's 
distribution). The results demonstrate the agreement with this hypothesis, though not 
large deviations are observed for small distances between levels. Such a tendency to the 
rise of some repulsion in regular region, apparently, is connected with a small admixture 
of chaotic component. 

The natural question about the role of numerical errors and, in particular, about 
the role of basis cutting off arises in careful studies of numerous examples confirming 
the hypothesis of the universal character of fluctuations of energy spectra. Even the 
special expression "the induced nonregularity" has appeared [06 |. The investigation of 



statistical properties of spectrum in three-levels model of Lipkin-Meshkov-Glick (LMG) 



[37J sheds a certain light on this problem. 

Let us consider the system of N particles on three levels, each of which is degener- 
ated with multiplicity iV . The Hamiltonian of this system in representation of second 
quantization has the following form 



H = Y, e* MT a+ a*, - -=V V (1 - S u ) V a+ or,, (36) 




where one-particle states (levels), marked by the indexes k (k — 0, 1, 2) , and n — 1, 2, ...N 
are different degenerate states belonging to one level,a + ,a; are fermionic creation and 
annihilation operators that obey the usual anticommutation relations. 

More traditional two-levels model LMG |38j has unique degree of freedom (it is the 
number of particles on top level) and consequently its classical analog does not assume 
chaotic behaviour. Introduction of an additional level is equivalent to the appearance of 
additional degree of freedom. Stationary states of the Hamiltonian fl3"6"|) can be marked 
by two indexes - the numbers of particles on the second and the third levels. As there is 
the finite number of ways of distribution of N particles on three levels, then the number 
of basis vectors in this model is finite. The problems of basis cutting off and estimation 
of arising errors are absent. 

Using coherent states 

\z) = exp z k a L a on |0> (37) 



v jfc=l n=l 
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where |0) is the state with all N particles on level with k = , we determine the classical 
Hamiltonian H c i as 

H d = ±-{z\H\z) =-l + ^ 1 2 (l- X ) + ^ 2 2 (2-x) + ^i(l + x) + ^(2 + x) + 



{Qi + (if - (pi +Pl) 2 ~ (Qi ~ Pi) (<?2 - P2) ~ 4qiq 2 piP2 

(38) 

Here x is the distance between levels, x = (N — l) v je , 

q k = Re— ^=,Pk = Im (39) 



1 + \zi\ + \z 2 \ y 1 + |^i I + \z 2 \ 

The equations of motion for q k and p k , obtained from variational principle are 
Hamiltonian. The classical analog of the Hamiltonian LMG is integrable at x = (V = 
0) . Changing the value of the constant of connection x we can observe the transition 
from regular motion to chaotic one. The spectral fluctuations for three-levels model 
LMG given in Fig. 22. These results are the serious argument in favour of the hypothesis 
of the universal character of fluctuations of energy spectra. 



1.5 R—C—R transition and the structure of QON wave function 

As we have seen in the previous section, the statistical properties of quantum spectrum 
of the Hamiltonian QON have been rigidly correlated with the type of classical motion. 
It is naturally to try to discover the analogous correlations in the structure of wave 
functions, i.e. to assume, that the form of the wave function for a semiclassical quantum 
state, associated with classical regular motion in the regions R\ or R 2 , is different 
from the one for chaotic region C . Notice that in contrast to the description of the 
energy spectra where we now have a sufficiently complete understanding of the spectral 
statistical properties we are still far from a corresponding complete knowledge of generic 
structure properties of the wave functions. Furthermore it should be pointed out that in 
analysis of QMCS at the level of energy spectra the principal role was given to statistical 
characteristics, i.e., the quantum chaos was treated as a property of the group states. 
The choice of a stationary wave function of the quantum system, which is chaotic in the 
classical limit, as an basic object of investigation relates the phenomenon of quantum 
chaos to an individual state. 

In contrast to spectrum the form of wave functions depends on the basis in what 
they are determined. Studying QMCS the three following representations are used more 
often 
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Figure 22: The spectral fluctuations in three-levels model LMG: a - distribution function 
of the nearest-neighbour spacings, b - the variance [81]. 
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1. The so called H -representation is: the representation of eigenf unctions {tp n } of 
integrable part H of total Hamiltonian H = H + V . The main objects of investigation 
in this case are the coefficients of expansion C mn of stationary functions ip m in basis {(p n } 
. Hq -representation is natural at realization of numerical calculations, as diagonalization 
of Hamiltonian H is realized more often just in this representation. 

2. Coordinate representation in which the behaviour of wave functions simply allows 
visual comparison with the picture of classical trajectories in coordinate space. 

3. Representation with the help of Wigner's functions [39] has a set of properties, 
common with classical function of distribution in phase space. 

As early as 1977 Berry [HO] assumed, that the form of the wave function ip for a semi- 
classical regular quantum state (associated with classical motion onaAf -dimensional 
torus in the 2N -dimensional phase space) was very different from the form of ip for an 
irregular state (associated with stochastic classical motion on all or part of the (2N — 1) 
-dimensional energy surface in phase space). For the regular wave functions the average 
probability density in the configuration space was seen to be determined by the projec- 
tion of the corresponding quantized invariant torus onto the configuration space, which 
implies the global order. The local structure is implied by the fact that the wave function 
is locally a superposition of a finite number of plane waves with the same wave number 
as determined by the classical momentum. In the opposite case for the chaotic wave 
functions the averaged (over small intervals of energy and coordinates) square of the 
eigenfunctions in the semiclassical limit h — > coincides with the projection of the clas- 
sical microcanonical distribution to the coordinate space. Its local structure is spanned 
by the superposition of infinitely many plane waves with random phases and equal wave 
numbers. The random phases might be justified by the classical ergodicity and this as- 
sumption immediately predicts locally the Gaussian randomness for the probability of 
amplitude distribution. Such structure of wave function is in well agreement with the 
picture of classical phase space: the classical trajectory homogeneously fill isoenergetic 
surface in condition corresponding to chaotic motion. By contrast, from the consid- 
eration of regular quantum state as analog of classical motion on torus, a conclusion 
should be done about the singularity (in limit h — > ) of wave function near caustics 
(boundaries of region of classical motion in coordinate space). 

Berry's hypothesis was subjected to the most complete test for billiards of different 



types and, in particular, for billiard stadium type [4T|. The amplitude of typical wave 



function of integrable circular billiard is negligible in classical prohibited region (from 
conservation of angular moment in circular billiard it should be that the arbitrary tra- 
jectory is enclosed between external and certain internal circle, the radius of which is 
determined by the angular moment) , whereas near the caustics it is maximal. Distri- 
bution |?/>| 2 for the case stadium, at which classical dynamics is stochastic, is strongly 
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distinguished from the integrable case. However, this distribution is not so homogeneous, 
as we could expect starting from the ergodicity of classical motion. 

Investigating QMCS in connection with properties of wave functions let us pass from 



billiards to Hamiltonian system of general position [42], f43|. For this purpose let us return 
to the object of our investigation: C$ v symmetric Hamiltonian QON. We begin to study 
the correlations of our interest by the topography of level curves of the stationary wave 
functions and, in particular, of nodal curves on which ipk{x, y) = . One of a nonrigorous 
criterion of stochasticity ||4| states that the system of nodal curves of regular wave 



function is a lattice of quasiortogonal curves or is similar to such a lattice. At the same 
time, the wave functions of chaotic states do not have such a representation. Fig. 23 
confirms that the structure of the lattice of nodal curves of separable wave functions 
undergoes a change in the R — C — R transition. The spatial structure of nodal curves 
for states from regions R\ and i?2 of regular classical motion is considerably simpler than 
this structure for states from region C of chaos. 

Correlations between the structure of wave function and the type of classical motion 
are also demonstrated in Fig. 24, in which the probability density \ip(x, y)\ 2 for states 
with numbers 111, 160 and 210 is represented. The squared module of the wave functions 
reproduces rather well transition from functions with clear internal structure (region Ri 
) to an irregular distribution (region C ) and the restoration of structure in the second 
regular region (i?2 ) • For the chosen technique, in which transition is traced for the wave 
function with fixed number (because of our use of the scaled Planck constant), a change 
in the wave function is associated only with R — C — R transition. 

Evolution of the wave functions in the process of R — C — R transition can be studied 
also in Ho -representation or more exactly in the representation of linear combinations of 
wave functions of two-dimensional harmonic oscillator with equal frequencies (see Section 
3.2) 

^ k = J2 C NLj\NLj) (40) 

N,L 

If one introduces the notion of distributivity of wave function in basis, then the 
criterion of stochasticity formulated by Nordholm and Rise [HH] even in 1974 states that 



with the degree of stochastization in the average arises the degree of distributivity of 
wave functions. It is clear that this criterion is a direct analog of the Berry's hypothesis 
for Hq - representation, if one interprets the number of basis state i = {NLj} as a 
discrete coordinate. Fig. 25 qualitatively confirm validity of this criterion. It can be 
seen from this figure that the states that correspond to regular motion are distributed 
in a relatively small number of basis states. At the same time, states corresponding to 
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k=111 k=16Q k ^ 210 




Figure 23: Nodal curves of the wave functions ipk (x,y) . The numbers k of the corre- 
sponding levels are shown over the figures. 
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Figure 24: Isolines of probability density \ipk(x, y)\ . The step between isolines is 
constant. The numbers k of the corresponding levels are shown over the figures. 
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Figure 25: Distribution of the coefficient C\ in the number i = {N,L,j} of the basis 
state. The numbers k of the corresponding levels are shown over the figures. 



chaotic motion (region C ) are distributed in a considerably larger number of basis state. 
In the latter case, the contributions from a large number of basis states in expansion 
p0| ) interfere; this results in a complex spatial structure of the wave function ipk {x, y) ■ 
For quantitative estimation of a degree of distributivity of wave functions, it is useful to 
introduce |46 , f47f the analog of usual thermodynamic entropy 

nk \ ^ I /^ik 

D ~ ~ / < \ ^NLj\ 



|2 in \r k I 2 

111 \ U NLj\ 



(41) 



Fig. 26 represents the dependence of entropy S from the number of state k for the 
first regular (Ri) and chaotic (C) regions. From this figure we can see, that the type 
of classical motion in large measure determines the character of this dependence: in 
chaotic region the entropy is practically constant, that corresponds to the high degree of 
distributivity, and in regular region the entropy is less in the value and has nonmonotone 
character. Nonmonotonicity of S is connected with rather high degree of localization of 
the wave function depending on the number of state. This effect and also the behaviour 
of the entropy in the second regular region R2 , will be studied in detail in the following 
section. 
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Figure 26: Entropy S as a function of the state number. The stright lines connecting 
the points correspond to shell classification according to N . 
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1.6 Evolution of shell structure in the process of R — C — R 
transition 



In sections 3.4 and 3.5 we have studied the change of statistical roperties of energy spectra 
and structure of wave functions in the process of R — C — R transition. However, the 
utility of the concept of the stochasticity extends further. Introduction of this concept 
in the nuclear theory made possible to take a fresh look JJJ at the old paradox ||: how 
one could reconcile the liquid drop - short mean free path - model of the nucleus with 
the independent particles - gas - like shell model. For the solution of this paradox within 
the limits of philosophy of simple chaos (see section 1.1 (part 1)) it is sufficient to assume 



1) When the nucleonic motion inside the nucleus is integrable, one expects to see 
strong shell effects in nuclear structure, quite well reproducible, for example, by the 
model of independent particles in the potential well. 

2) When nucleonic motion is chaotic, one expects smooth, statistical, Thomas-Fermi, 
Droplet model approaches to be good approximations. This is because phase space is 
much more uniformly filled with the Boltzmann fog in this case. At such approach the 
elucidation of the mechanism of destruction of shell effects in the process of the transition 
regularity- chaos plays the key role. More appropriately the problem can be formulated 
in the following way J?|: how do shell dissolve with deviations from integrability or, 
conversely, how do incipient shell effects emerge as the system first begins to feel its 
proximity to an integrable situation? 

As it has been mentioned many times, the finite motion of integrable Hamiltonian 
system with N -degrees of freedom is, in general, conditionally periodic, and the phase 
trajectories lie on N -dimensional tori. In the variables action-angle (J, 9) the Hamil- 
tonian is cyclic with respect to angle variables, H = Hq (I) . Even at the beginning 
of the century, Poincare called the main problem of dynamics the problem about the 
perturbation of conditionally periodic motion in the system defined by the Hamiltonian 



where e is a small parameter. The essential step in the solution of this problem was the 
KAM theorem, asserting, that at switching on nonintegrable perturbation, the majority 
of nonresonance tori, i.e. the tori for which 



is conserved, distinguishing from the nonperturbed cases only by small (to the extent 
of smallness of e ) deformation. As we have seen in section 2.3(part 1), at definite 








(42) 




(43) 
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conditions the KAM formalism allows to remove out of Hamiltonian the members de- 
pending on angle, using the convergent sequence of canonical transformations 0. When 
it is succeed, we immediately find that perturbed motion lies on rather deformed tori so 
that trajectories, generated by perturbed Hamiltonian, remain quasiperiodic. In other 
words the KAM theorem reflects an important peculiarity of classical integrable sys- 
tems to conserve regular behaviour even at rather strong nonintegrable perturbation. In 
the problem of our interest, concerning the destruction of shell structure of quantum 
spectrum, the KAM theorem also will be able to play an important role. Considering 
the residual nucleon-nucleon interaction as nonintegrable addition to selfconsistent field, 
obtained, for example, in Hartree-Fock approximation, one can try to connect the de- 
struction of shell with the deviation from integrability. The existence of shell structure 
at rather strong residual interaction (or at large deformation) can be obligated to rigid- 
ity of KAM tori, contributing to survival of regular behaviour. Such assumption seems 
rather natural, especially if one takes into account, that the procedure of quasiclassical 
quantization itself |48], [49] is indirectly based on convergence of the same sequence of 
canonical transformations, as the KAM theorem. 

The aim of this section is to trace the evolution of the structure of QON Hamiltonian. 
In numerical calculations of this section it will be more convenient for us to use nonscale 
version of the Hamiltonian (16.1). By nonperturbed Hamiltonian Ho we will mean 
the Hamiltonian of two-dimensional harmonic oscillator with equal frequencies: H$ = 
H (a = 1, b = 0, c = 0) . Its degenerate equidistant spectrum is well known. At switching 
on perturbation, the degeneration disappears and the shell structure forms. The number 
of states, for example, the states of E - type (the numerical results, represented below, 
are relative to the states precisely of this type, though analogous results take place for 
the states belonging to another nonreducible representations of C$ v - group) for the given 
quantum number N is equal to \ (N 1 + Mod (N, 2)) , where = | (2iV + Mod (N, 3)) 
, and Mod (N,M) is a remainder from division of N by M . 

It is obvious, that eigenfunctions of exact Hamiltonian QON are not the eigenfunc- 
tions of operators N and L already. Nevertheless, as numerical calculations show even 
at rather large nonlinearity, one can use the quantum number N and L for the classifi- 
cation of wave functions. The measure of nonlinearity, at which such classification (i.e. 
the existence of shell structure) remains reasonable, is connected with quasicrossing of 
neighbouring levels. Under the quasicrossing we shall understand the approach of levels 
up to the distances of degree of precision of numerical calculations (the detailed analysis 
of quasicrossing see lower). 

The dependence of energy spectra of Hamiltonian QON from the parameter b for the 
values W = 3.9 and W = 13 are represented in Fig.27a,b. As it is seen from Fig. 27a. for 
the PES with W = 13 at the approaching to the line of critical energy of the transition 
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Figure 27: a. Energy spectra of Hamiltonian (16.1) depending on values of parameter 
b for region C (W = 13) . Points mark quasicrossinds. Continuous line - dependence 
E cr from b . Dashed line shows the beginning of the region of quasicrossings. Arrow 
shows point of quasicrossing of the levels with k = 40 and 41 . b. Energy spectra of 
Hamiltonian (16.1) depending on values of parameter b for region Ri (W = 3.9) . There 
are no quasicrossing of the levels. 



to chaos, defined according to the criterion of negative curvature, the destruction of 
shell structure, which we understand as multiple beginning of quasicrossing, takes place. 
At the same time, for the PES with W = 3.9 in Fig.27b (for which as we have shown 
in Section 2. 4 (part 1), the local instability is absent and the motion is regular at all 
energies) the quasicrossings are absent even at larger nonlinearity than in Fig. 27a for 
W = 13 . 

The destruction of shell structure can be traced for the wave functions, using intro- 
duced in previous section, the analog of thermodynamic entropy S . The character 
of changes of entropy in regions R± and R2 , corresponding to regular classical motion, 
correlates with the transition from shell to shell (see Fig. 28). Two effects are observed 
in the region C , corresponding to chaotic classical motion. Firstly, quasiperiodical 
dependence of entropy from energy is violated, that testifies to destruction of sell struc- 
ture. Secondly, the monotone growth of entropy is observed on average with going out 
on plateau corresponding to the entropy of random sequence at the energies essentially 
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Figure 28: Entropy S and quantum numbers L,N as the functions of state number k. 



exceeding critical energy. 

Average values of operator N , given in Fig. 29, and calculated on stationary wave 
functions of exact Hamiltonian, describe the dynamics of change of the shell structure. 
From this Figure we can see, that the minimum deviation of (N) from N is observed in 
regions R\ and R2 , while this deviation is considerably greater in the stochastic region 
C ; moreover in the average it monotonically increases with increase of energy. It is 
obvious, that at energies, for which N — (N) ^ 1 , the destruction of quantum analogs 
of classical integrals leads to the fact that the classification of levels with the help of 
quantum numbers N and L loses sense. 

Notice, that we could observe restoration of shell structure at high energies in the 
process of R — C transition only due to the optimization of basis frequency (see Section 
1.2), that is equivalent to the possibility of diagonalization of matrices of considerably 
higher dimension. 

Now let us return to the analysis of quasicrossings of levels arising near critical energy 



of the transition to chaos. As it is known [50|, for the set of Hamiltonians depending 



33 



A'. <iV> 



20- 



10- 




20 



10 



t 



■v" J.' 



20- 



10 



40 



SO 



Figure 29: Dependence of N (lines) and (JV) points) of state numbers k 
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on two parameters A = (b, c) and that is invariant with respect to inversion of time 
in the neighbourhood of degeneration point A* = (b*,c*) , the difference of two terms 
AE — E\ — E 2 is determined by the expression 

AE (b, c) = [A(b- b*f + B{b- b*) (c -c*) + C{c- c*) 2 ] 11 (44) 

where coefficients A, B, C are the functions of components Va-^ in the point A* . The 
characteristics of nuclear shape [H8| in the problems of nuclear spectroscopy more often 



play the role of parameters. The wave functions of intersecting levels near degeneration 
point can be represented in the following form 

\Ex (A)) = cosxi |£i (A*)) + sinxi \E 2 (A*)) 

\E 2 {X)) = cosx 2 |£i(A*)> + sinx2|£ 2 (A*)> 1 j 

where angles of mixing xi,2 are determined by the relation 

e**,, = ± (TX)+iiSX) 

[(TA) 2 + (5A) 2 ] 1/2 ' (46) 
(TA) « H u (A) - H 22 (A) , (SX) « 2H 21 (A) 

It is easy to see, that changing A — > —A the "exchange" of wave functions of approaching 
levels takes place \E\ (A*)) <-> \E 2 (A*)) . Such test can be effectively used for the analysis 
of nature of quasicrossing points. 

We trace, in particular, the variations of functions with K = 40 and K = 41 , which 
undergo quasicrossing in the neig hbourhood b* w 9.8 • 10~ 2 , c* w 7.4 ■ 10" 4 (recall, 
that under this term we understand the approaching of levels up to distances of an order 
of precision of numerical calculations). The sections of ty(x = 0,y) and coefficients of 
expansions C% L j of wave functions of considering states before and after quasicrossing 
point are represented in Figs 30 and 31. The similar characteristics of the states with 
K = 39 and K = 42 which do not undergo quasicrossing at this values of parameters 
are represented for comparison at the same figure. 

As it is seen in represented figures the wave functions are not practically changed for 
the states 39 and 42, at the same time the exchange of wave functions is observed for 
the states 40 and 41 in correspondence with test described above. Analogous situation 
takes place for other quasicrossings . 
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Figure 30: Wave functions ip(x = 0,y) for the states with k = 39,40,41,42 in the 
neighbourhood of quasicrossing of the states with k = 40 and 41 . Left b < b* , right 
b > b* (point of quasicrossing: b* = 0.098; c* = 0.00074 ). 



1.7 Wave packet dynamics 

The investigation of the time evolution of nonstationary states, i.e. of packets in quan- 
tum systems classical analogs of which admit chaotic behaviour, gives an important 
information about QMCS. The localized quantum wave packet (QWP) is the closest 
analog of the point in phase space, which describes the state of classical system. How- 
ever, the correspondence between localized QWP and classical point particle in chaotic 
region is broken down in a very short time interval. Let us explain it by using Taka- 
hashi arguments [52]. We take two localized QWP (x) and ^2 {%) which are put in 
the chaotic region at initial time to be slightly different from each other so that the 
difference between and (\&2| £ 1^2) ( or an d (^l P ^2) ) is very 

small. We assume that in chaotic region the localized QWP does not either extend in a 

certain time interval of the order J% like that in regular region. Thence, following the 
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Figure 31: Same as Figure 30, but for the coefficients of the expansion of wave functions 



l*(0| 2 = 



, i , - number of basis state. 



Ehrenfest theorem, the packets and ^2 move as classical particles and the distance 
between £ and (^l x [^2) (or and (^l P l^) ) is increasing expo- 

nentially in time. Let us consider the superposition ^ = ^1 + ^2 which also becomes 
a localized QWP in initial state. From the assumption, we can expect that QWP ^ 
does not extend in a certain time interval of the order ^ I yf^ . However, considering the 

exponential increment of the distance between x \^\) and (^ 2 | x |\l/ 2 ) , ^ extends 
exponentially in chaotic region and does not behave as a classical particle. This result 
is inconsistent with initial assumption and implies that in chaotic region the localized 
QWPs (i.e. ^1 , ^>2 and \1/ ) extend exponentially like a classical probability distribution 
in the first stage of time development. 

In order to describe such unusual behaviour of QWP one should address to the con- 
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cept of a quantum-mechanical phase space. There are a few well-established schemes to 



introduce phase-space variables in quantum mechanics [39, 53, 64, 55]. In the present 



study we shall follow the procedure proposed by Weissman and Jortner |pb[ . Let us con- 
sider an initially localized QWP \I/ characterized by the coordinates q and the momenta 

P , 



q = (*| q |*) 
p = (^\p\^) 



(47) 



Now we introduce the coherent states \p, q) |5 
are given by Gaussian wave packets 



which in the coordinate x-representation, 



N 



(x\p,q) = Y[ (tto-J) /4 exp 



\ x j ( h 1 
2a) 



+ 



iPjXj ip-jq-j 



h 



2h 



(48) 



In the study of a system of iV coupled harmonic oscillators, it is convenient to choose 
for constants Oj the rms zero-point displacements 



cr-i 



y 1 



V2 



(49) 



where mj are masses and uuj are frequencies of the uncoupled oscillators. With this 
choice of the aj , the coherent states \a) = \p, q) become the eigenstates of the harmonic 
oscillator annihilation operators aj 



where a complex variable a 



Q'; 



cij | a) = otj | a) 



1 f Qj , -°i 



(50) 



(51) 



Using these coherent states, it is possible to introduce the following quantum - me- 
chanical phase-space density 



(52) 



where \1/ is any general wave packet. This quantum- mechanical, coherent-state, phase- 
space density may be regarded as a quantum analog of the classical phase-space density, 
since it satisfies an equation of motion whose leading term (when expanded in powers 
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of H ), corresponds to the classical Liuville equation [p^j . The stationary phase space 
densities (p, g) are the squares of the projections of the eigenstates on the coherent 
state, 

p E (p,q) = \(a\E)\ 2 (53) 

where |a) is given in terms of Gaussian wave packet (f48f ) while the eigenstate \E) is 
given by Eq.(|J). Using the well-known expressions for scalar products (a \ NL) |54|] we 
finally obtain 

Pe (p, q) = - exp 

E^NL I *nx *n— , ■ *n— *n+\ 
— [o. + a_ + j a + a_ J 

n,l (n_\ n + \) 7 2 

where a± = (a 2 =F with aj )2 = ^ (gi, 2 + zpi, 2 ) , n + = , n_ = • The 
phase-space density p# (p, g) is a function of the four real variables pi , gi and p 2 , g 2 • 
We can get the contour maps of pE (pi, gi,P2, 92) in the g 2 x p 2 plane, taking gi = and 
calculating pi from the relation 



~\ (l«+| 2 + I"-! 2 ) 



ff(gi = 0,pi,g 2> p 2 )=E (55) 

Quantum Poincare maps (QPM) |56| , obtained by such way, constitute the quantum 
analogs of the classical Poincare maps and can be used for the search of QMCS both in 
the structure of wave functions of stationary states and in the dynamics of wave packets. 

Now we shall consider the time evolution of a wave packet which is initially in a 
coherent state 



|*(t = 0)) = |or) (56) 
The time evolution of such an initially coherent wave packet is given by 

\^(t)) = J2\Ek)(E k \a)e- lEkt (57) 

k 

The survival probabilityp (t) of finding the system in its initial state is 

p(t) = \g a (t)\ 2 (58) 



39 



where g a (t) is the overlap of \I/ (t) with the initial state 

g a (t) = (a | * (t)> = X) I I »> I' ( 59 ) 



Utilizing the definition of the stationary phase-space density (p3|), we can write 

9a {t)=Y,PE k e- lEkt (60) 



Equation ( |6"0"D implies that dynamics is determinated by the spectrum of the initial 
coherent state |a) . 

Weissman and Jortner |5?J] have observed for Henon-Heiles Hamiltonian (QON Hamil- 
tonian with m = l;a=l;c = 0) two limiting types of QWP dynamics of initially 
coherent Gaussian wave packets, which correspond to quasiperiodic time evolution and 
to rapid decay of the initial state population probability Quasiperiodic time evolution 
is exhibited by wave packets initially located in regular region, while rapid decay of the 
initial state population probability is revealed by those wave packets that are initially 
placed in irregular regions. 

Ben-Tal and Moiseyev p5| calculated the survival probability p (t) for an initial com- 
plex Gaussian wave packets by the Lanczos recursion method. This method has a prac- 
tical value since it requires nN 2 numerical operations (n is number of Lanczos recursion, 
N is dimension of the original Hamiltonian matrix) rather than N 3 required to calculate 
the eigenvectors H . For bounded time calculations n is much smaller than N . 

Let us turn to the consideration of the dynamics of QWP in the PES with some 
local minima (W > 16) . Transitions between different local minima can be divided 
into induced (the excitation energy exceeds the value of potential barrier) and tunnel 
transitions. The latter ones are subdivided into transitions from discrete spectrum into 
continuous spectrum (for example, a -decay, spontaneous division) and from discrete 
spectrum into discrete spectrum (for example, transitions between isomeric states). The 
process of tunneling across a multidimensional potential barrier, when initial and final 
states are in discrete spectrum, is the most complicated and almost noninvestigated 
problem. 

The time evolution of wave packet is more often studied by two methods [57]: ei- 
ther by direct numerical integration of the Schrodinger time-dependent equation with 
corresponding initial condition \& (r, t — 0) , or by expansion of the packet ^ (r, t) in 
the eigenfunctions of the stationary problem. The first method has some lacks, e.g., 
the complication of interpretation of the obtained results and the necessity of separating 
the contributions from subbarrier and tunnel transitions for the packet of an arbitrary 
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shape. These difficulties can be avoided providing that the subbarrier part of the spec- 
trum E n (E n < Uq, Uq -the height of barrier) and the corresponding stationary wave 
functions ip n (r) are known. A pure tunnel dynamics will take place for the packets 
representable in the form 

(f, t) = J2 C n e-i Ent i> n (r) ,E n <U (61) 

n 

where 

C n = J * (r, t = 0) i) n (f) df (62) 

The probability p R (t) of finding the particle at the moment of time t in certain local 
minimum R is 

p R (t) = J \^(r-t)\ 2 df=Y,C* m C n e 11 ^^ I r m (r)^P n (^df (63) 

m,n 

or in the two level approximation 

P R (*) = P R (0) - 4CiC 2 sin 2 / ^i-^m /• ^ ( f) ^ ( f) rff (g4) 



Let us introduce |^8[ the value 

p fi = max/ (t) (65) 



vt 



which is the maximum probability of finding the particle in the certain local minimum 
R , if initially it was localized in the arbitrary minimum. If the number of local minima 
is more than two, then of independent interest is 



p Ro = mmp Ro (t) (66) 

Vi 

i.e., minimum probability of finding the wave packet in the minimum R corresponding 
to its initial localization. 

Intuitively, we may suggest that p R m 1 if the initial minimum is local, and the 
final one is absolute. However, the results |5£| obtained for the simplest one-dimensional 
models (asymmetric double wells of different shapes) are inconsistent with the intuitive 
expectations. The probability of tunneling from the local minimum to the absolute one 
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Figure 32: Dependence of subbarier energy levels Ei for double asymmetric one- 
dimensional rectangular well with infinite external walls (a) and p R , (b) from d . Width 
of well is equal to 3 , width of barrier - 1 , height of barrier - 2 . 



N(I! 
1000 



500 







0,5 0,6 0,7 0,8 0,9 1,0 I 



Figure 33: Distribution of subbarier wave functions iV (J) 
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resonantly depends on the potential parameters. Fig. 32 gives the dependence of p R from 
well depth displacement d. It is seen that at an arbitrary asymmetry p R ^ 1 . 

The resonant behaviour of p R becomes clear if one considers the spatial structure 
of the subbarrier wave functions. For a sufficiently wide barrier at an arbitrary asym- 
metry, the subbarrier wave functions are largely localized in separate minima. Fig. 33 
gives the distribution iV (/) of the subbarrier wave functions as their degree of localiza- 
tion I n = J \ijj n (f)\ 2 df for an accidental choice of the well depth displacement d. The 

R 

derealization takes place only in the vicinity of the quasicrossing of level. The degree 
of this derealization is directly depends on the distance between the interacting levels. 
Obviously, the QWP, in which the components localized in the certain minimum are 
dominating, cannot effectively tunnel to the neighbouring minimum. This just explains 
the stringent correlation between the p R minima and the level quasicrossing (see Fig. 
32.) 

Now the question arises if a similar correlation between the level quasicrossings, the 
derealization of wave functions and resonant tunneling persists in the two-dimensional 
case. To give an answer to this question let us turn to the results of numerical investi- 
gation of stationary wave functions of Hamiltonian QON in region W > 16 |59]. Recall 
that three {C% v - symmetry) identical additional minima appear at W > 16 apart from 
the central minimum. The central minimum exceeds the lateral ones in depth in the 
region 16 < W < 18. At W > 18 , the central minimum becomes the local one. In this 
region of parameters the procedure of diagonalization of Hamiltonian QON in oscillator 
basis essentially complicates. It is connected with the fact, that the basis of Hamilto- 
nian, the potential of which has the unique minimum, is used for the diagonalization of 
Hamiltonian possessing complex topology of the PES. 

In addition to large dimension of the basis it is necessary that the basis wave functions 
have sufficient value in the region of lateral minima. It is achieved by the optimization 
of oscillator frequency of basis uq . For the values of parameters used in the calcula- 
tions (W = 17.8; b = 0.17) the value of oscillator frequency loq = 0.2 . Fig. 34 gives low 
eigenvalues of E and A 1 types depending on the basis dimension. We can see, that it is 
possible to get the saturation in basis at the dimension of submatrices N ~ 10 3 even for 
low states localized in the lateral minima (dashed lines). 

We can also use introduced above the analog of thermodynamic entropy Sk for esti- 
mation of degree of saturation in basis. Fig. 35 gives the Sk values for the states of Ai 
type for the dimension of submatrices 408 and 690. Increasing of the basis does not lead 
to sufficient changes of values Sf. for the states with energy up to saddle energy E s ■ 

The states, localized in the central or in the lateral minima, have essentially different 
distributivity of coefficients C\ {i = NLj}(see Fig. 36a,b) and thus different entropies: 
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Figure 34: Dependence of energy spectrum E -type (a) at b = 
(b) at b = 0.17 from the dimension of submatrices iV for W 
correspond to the states localized in peripheral minima. 



0.173798 and A x -type 
= 17.8 . Dashed lines 




Figure 35: Dependence of Sk from number of eigenstate k {A\ - type) for the dimension 
of submatrices 408 (dark point) and 690 (crosses) at W = 17.8 and b = 0.17 ; E s 
corresponds to saddle energy and E cr - the classical critical energy of the transition from 
regular type of motion to chaotic one. 
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Figure 36: Distribution of coefficients |C;| by number of basis state i = {N,L,j} (E - 
type) for the localized states with k = 3 (a) and k = 4 (b) at b = (U7379 (W = 17.8) 
and for the delocalized states with k — 3 (c) and k — 4 (d) in the point of quasicrossing 
6 = 0.1737924 (W = 17.8) . 
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Figure 37: Dependence of energy levels Ei of Hamiltonian (16.1) from the parameter b 
for W = 17.8 : a - spectrum of E - type, b - spectrum of Ai - type. 

states, localized in the central minimum, have less entropy). In the neighbourhood of the 
points of level quasicrossings, the derealization of wave functions, corresponding to these 
levels, takes place; these wave functions possess close distributivity of coefficients C\ (see 
Fig. 36.c.d). Fig. 37 shows the subbarrier part of the energy spectrum obtained by the 
diagonalization. As it is easy to see the tunneling of the wave packet, composed of the 
subbarrier wave functions, can be described in the two-level approximation. Indeed, there 
are approximately 10 level quasicrossings of A\ and E - types, where the nonlinearity 
parameter changes are of the order of 10~ 2 . The effective half-width of the overlap 
integral in ( |5~4"[ ) is about 10~ 5 (see Fig. 38). Hence, all nondiagonal elements will be 
close to zero with the overwhelming probability in the matrix a mn = J tfrni* (r) ip n (r) df 

R 

for an arbitrary nonlinearity parameter (e.g., b ). Two appreciable different from zero 
nondiagonal matrix elements (two-level approximation) appear only in the vicinity of 
quasicrossings. The probability of double quasicrossings at a fixed nonlinearity parameter 
is almost excluded. This probability is by two or three orders lower than that of rather 
rare ( 10~ 3 ) single quasicrossings. 

So, now we can give the answer to the question we have posed above. The stringiest 
correlations between level quasicrossings, derealization of wave functions and resonant 
tunneling across the potential barrier take place in the two-dimensional case (and, most 
likely, in the multi-dimensional one too). 

The existence of the mixed state for the many-well potentials must be well manifested 
itself in dynamics of the QWP. Preexponential factor of the tunnel amplitude depends on 
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Figure 38: a - quasicrossing of energy levels of E - type with k — 3, 4 for Hamiltonian 
(16.1); b - localization of wave functions of states with k = 3,4 in the central well at 
different b ; c - dependence p R from b . 



the type of classical motion and consequently we expect to observe the asymmetry of the 
effective barrier penetration in the mixed. This purely classical effect may be observed 
only if the uncertainty in level energy is comparable with the average distance between 
levels and the system does not yet "feels" that the spectrum is discrete. This determines 
the time scale for which observation of effect is possible. It is the same time scale on 
which transition from classical linear diffusive increasing of energy to quasiperiodical 



quantum evolution is observed [60 



2 Summary and open questions 
2.1 Summary 

In this review we have presented a complete description of classical dynamics gener- 
ated by the Hamiltonian of quadrupole oscillations and identified those peculiarities of 
quantum dynamics which can be interpret as QMCS. 

We have shown an intimate connection between dynamics features and geometry 
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of the PES. Interpretation of negative curvature of the PES as the source of the local 
instability allowed to correctly predict the critical energy of the transition to chaos for 
one-well potentials. 

Particular attention has been given to investigation of classical dynamics in the pa- 
rameter region according to the PES with a few local minima. As is shown, one of the 
main peculiarities of the many-well Hamiltonians is the mixed state: at one at the same 
energy in a different local minima diverse dynamical regimes (regular or chaotic) are 
realized. Until the present time the most theoretical and numerical works have been 
focused on the behaviour of so-called billiards or Hamiltonian systems with the simplest 
topology of the PES. Evidently it is not enough for understanding of the real many-body 
systems with the complicated PES for which the mixed state is a situation of general 
position. Therefore the present review should be read as one of the indispensable steps 
on the way to transition from description of the model systems to direct consideration 
of much more realistic systems. 

A new methodology for determination of the critical energy of the transition to chaos 
based on the investigation of correlations between dynamic features and geometry of the 
PES is developed for many-well potentials. It is more simpler than criteria of transition 
to chaos connected with one or other version of overlap resonance's criterion. At the 
energies below the critical one approximate integrals of motion for every of local minima 
are derived by the Birkhoff-Gustavson normal form technique. 

We have numerically demonstrated that for potential with a localized unstable region 
(in particular localized region of negative Gaussian curvature) regular motion restores at 
high energy, i.e. transition regularity-chaos-regularity takes place for these potentials. 

QMCS in dynamics of QON have been the center of our consideration. The variations 
of statistical properties of energy spectrum in the process of R — C — R transition have 
been studied in detail. For the chaotic region all the analyzed statistical characteristics 
are seen to be in good agreement with the GOE predictions. For regular regions our 
calculations are consistent with the hypothesis of the universal character of energy spec- 
trum fluctuations. Though for small level spacings we observed some deviations which 
were probably due to a small admixture of the chaotic component. 

We proved that the type of classical motion is correlated with the structure of the 
stationary wave functions of highly excited states in the R—C—R transition. Correlations 
were found both in the coordinate space (the lattice of nodal curves and the distribution 
of the probability density) and in the Hilbert space associated with the integrable part 
of Hamiltonian (the distribution of the wave functions in the oscillator basis and the 
entropy of individual eigenstates). Calculations with the scaled Planck constant, which 
make it possible to obtain wave functions with equal quantum numbers and energies that 
correspond to different types of classical motion, enabled us to unambiguously isolate 
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correlation effects in the structure of wave functions. 

The Hamiltonian of QON was used as an example to study the shell structure de- 
struction induced by the increase of nonintegrable perturbation term which models resid- 
ual nucleon-nucleon interaction. In the vicinity of the classical critical energy value it 
was observed multiple formation of quasicrossings of the energy levels, violation of the 
quasiperiodical energy dependence of the entropy, and the increase of the average value 
fluctuations of the operators used to classify the eigenstate of the integrable problem. 

The optimization of the basis frequency and, as a consequence, the possibility of 
diagonalization matrices of higher dimension made it possible to trace the restoration of 
the shell structure in the transition C — R2 from chaos to regularity. 



2.2 Open questions 

The problem of studying of QMCS in Hamiltonian systems with a few local minima 
represents almost a noninvestigated region. We must take into account that, apart from 
conceptual difficulties, construction of spectrum and eigenfunctions in this case is a very 
difficult computational problem. Calculation of energy splitting DE due to tunneling 
is one of the oldest problem of quantum mechanics. One dimensional problem for the 
most part is understood. A complication arises in the multidimensional problem, which 
might not have been expected until recently. This is because the semiclassical wave 
functions which determine energy splitting are sensitive to the nature of the classical 
motion. Wilkinson [104] calculated the energy splitting due to tunneling between a pair 
of quantum states which correspond to classical motion on tori in phase space. The 
case where both wave functions correspond to classical chaotic motion were discussed by 
Wilkinson and Hannay [105]. Problem of the wave function structure for the mixed state 
remains to be solved. The configuration space in this situation breaks up naturally into 
different regions (separate local minima), in each of which different dynamics regimes are 
realized. Thus, there is a need of a calculating scheme, which will enable us to use the 
partial information about each isolated region to obtain solution of the full problem. We 
believe that the path decomposition expansion of Auerbach and Kivelson [B3] will prove 
to be useful for solving 2 — d tunneling problem in the mixed state. This formalism allows 
to express the full time evolution operator as a time convolution and surface integrations 
of products of restricted Green functions, each of which involves the sum over paths that 
are limited to different regions of configuration space. Even for complicated nonseparable 
potential the qualitative behaviour is readily inferred and quantitative solutions can be 
obtained from knowledge of the classical dynamics. 

In conclusion, we note that an important problem concerning the role of periodic orbit 
in the structure of wave functions was not considered in this study. In the semiclassical 
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limit quantum state must be determined by the invariant sets formed by classical tra- 
jectories. These sets cover the total energy surface for systems that are chaotic in the 
classical limit. The approach in which the energy surface as a whole plays a dominant 
role is confirmed by the chaotic structure of nodal curves and by the random distribution 
of the probability density of the wave functions of highly exited states. However, the 
recent findings of Heller |B3], who showed that wave functions can have the so-called 



scars corresponding to a high concentration of probability density near unstable periodic 
orbits, demonstrate that such an approach is not complete. The reason for this is that, 
although the measure of periodic trajectories is zero, their contribution is essentially sin- 
gular, in contrast to the smooth contribution from the energy surface. These problems 
require special analysis for the Hamiltonian under study. 
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